A Significant Problem With Using the Amati Relation for 

Cosmological Purposes 

Andrew C. Collazzi, Bradley E. Schaefer 
Dept. of Physics and Astronomy, Louisiana State University, Baton Rouge, LA 70803 

Adam Goldstein, Robert D. Preece 
Dept. of Physics, University of Alabama, Huntsville, AL 35899 

ABSTRACT 

We consider the distribution of many samples of Gamma-Ray Bursts (GRBs) 
when plotted in a diagram with their bolometric fluence (Sboio) versus the ob- 
served photon energy of peak spectral flux {Epeak,obs)- In this diagram, all bursts 
that obey the Amati relation (a luminosity relation where the total burst energy 
has a power law relation to Epeak,obs) must lie above some limiting line, although 
observational scatter is expected to be substantial. We confirm that early bursts 
with spectroscopic redshifts are consistent with this Amati limit. But, we find 
that the bursts from BATSE, Swift, Suzaku, and Konus are all greatly in viola- 
tion of the Amati limit, and this is true whether or not the bursts have measured 
spectroscopic redshifts. That is, the Amati relation has definitely failed. In the 
Shoio — Epeak,obs diagram, we find that every satellite has a greatly different dis- 
tribution. This requires that selection effects are dominating these distributions, 
which we quantitatively identify. For detector selections, the trigger threshold 
and the threshold for the burst to get a measured Epeak,obs combine to make a 
diagonal cutoff with the position of this cutoff varying greatly detector to detec- 
tor. For selection effects due to the intrinsic properties of the burst population, 
the distribution of -Epeafc.ofes makes for bursts with low and high values to be rare, 
while the fluence distribution makes for bright bursts being relatively uncommon. 
For a detector with a high threshold, the combination of these selection effects 
serves to allow only bursts within a region along the Amati limit line to be mea- 
sured, and these bursts will then appear to follow an Amati relation. Therefore, 
the Amati relation is an artifact of selection effects within the burst population 
and the detector. As such, the Amati relation should not be used for cosmologi- 
cal tasks. This failure of the Amati relation is in no way prejudicial against the 
other luminosity relations. 

Subject headings: gamma-ray burst: general, gamma-ray: stars 
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1. Background 

Luminosity relations are a tool for which we can connect a measurable quantity (called 
a luminosity indicator) from a long-duration Gamma-Ray Burst (GRB) to the burst's peak 
luminosity or total energy. The inverse square law and the observed burst brightness can 
then be used to determine the distance to the burst. In doing this, GRBs become a 'standard- 
candle' (in the same sense as Cepheid variables and Type la supernovae) where observed 
properties can be used to determine luminosity and then distance. GRBs are seen out to very 
high distances, at least to redshift z=8.2 (Tanvir ct al. 2009), and so the GRB luminosity 
relations offer a unique and powerful opportunity for new cosmology. 

Currently, GRBs have seven well-estabhshed luminosity relations (Fenimore & Ramirez- 
Ruiz 2000; Norris et al. 2000; Amati et al. 2002; Schaefer 2002; Schaefer 2003; Ghirlanda 
et al. 2004). The luminosity indicators are the spectral lag time {riag, Norris et al. 2000), 
'variability' (V, Fenimore & Ramirez-Ruiz 2000), the peak of the uF^, power spectrum {Epgak', 
Schaefer 2003), the minimum rise time in the light curve [t^it] Schaefer 2002), and the 
number of peaks in the light curve {Npeak, Schaefer 2002). Five of the luminosity relations 
connect the luminosity indicators to the burst's peak luminosity. Two luminosity relations 
instead connect to the total energy of the burst, with the physics relating to the burst 
energetics instead of the conditions at the time of peak luminosity. These are: the (A) 
Amati relation (Amati et al. 2002;2006), which connects -Bpeafe to the total energy of the 
burst assuming for isotropic emission {E^^iso), and (B) the Ghirlanda relation (Ghirlanda et 
al. 2004) which connects Ep^ak to the total energy of the burst, but with a correction for the 
beaming factor (-Ffeeam), resulting in the energy emission {E^). In addition, there are a variety 
of newer luminosity relations have been proposed, but not yet fully tested and confirmed by 
the community (e.g. Dainotti et al. 2008; 2010; 2011) and while others have proven to not 
be viable. An example of this is the Firmani relation (Firmani et al. 2006), which was found 
to be no significant improvement over existing relations (CoUazzi & Schaefer 2008). Physical 
explanations for a variety of the established relations have been made (e.g. Kobayashi et 
al. 2002; Meszaros et al. 2002; Schaefer 2003; Eichler & Levinson 2004; Liang et al. 2004; 
Schaefer 2004; Rees & Meszaros 2005; Giannios & Spruit 2007; Thompson et al. 2007). 

Recently, CoUazzi et al. (2011) performed an exhaustive study on the sources of error 
on Epgak- This study showed a considerable amount of scatter was hidden in how Ep^ak is 
measured, which was much larger than that from the reported Poisson errors alone. The 
sources of this scatter included: (A) the choices of different analysts, (B) which Epeak is 
measured (e.g. the time-integrated or time-resolve Epeak), and (C) the detector response 
matrix, all in addition to the regular Poisson statistical error. This scatter can be as large 
as 0.43 in log space and has a typical value of 0.24. This scatter can explain the scatter seen 



- 3 - 



in the luminosity relations that use Ep^ak- 

Still, even the currently accepted luminosity relations have their drawbacks. The best 
(i.e. the tightest) of these relations, the Ghirlanda relation, can only be applied if there is 
an observed jet break. Jet breaks are a well understood phenomena (Rhoads 1997; Sari et 
al. 1999). Measuring a jet break is fairly difficult for a variety of reasons, and has only been 
observed in a small percentage of bursts. Melandri et al. (2008) and Kocevski & Butler 
(2008) have pointed out problems in identifying these jet breaks with the X-ray data. 

Most notably, the Amati relation has been criticized for several reasons. Li (2006) orig- 
inated a test that demonstrated an ambiguity when the measured properties are used to 
determine the redshift. This result was confirmed by other groups later (e.g. Schaefer & 
Collazzi 2007). This ambiguity also exists for the Ghirlanda relation, but at a substantially 
higher redshift. This ambiguity does not exist for any of the other confirmed luminosity rela- 
tions, and so the ambiguity problem can go away if multiple relations are used to determine 
redshift. This problem does not affect current work on the ORB Hubble Diagram because 
those luminosity indicators were derived from a known redshift obtained spectroscopically. 

Another major criticism towards the Amati relation came from Nakar and Piran (2005). 
In this work, Nakar and Piran developed a test specifically for the Amati relation, the beauty 
of the test being that a redshift was not needed. This test has since been generalized in several 
independent investigations (e.g. Band & Preece 2005; Schaefer & Collazzi 2007; Goldstein 
et al. 2010). They combined the Amati relation (equation [T]) with the inverse square law for 
fiuences (equation |2]) to eliminate E^^iso (equation [3]) . 

E.^iso = (9.2 X 10^^ erg keV'^-o^) [E^ (i + ^)]2-04 (i) 
_ indlSboio , . 

Sbolo (1 + ^)=' °^(9.2 X 1047 erg keV-^-o^) ^ ^ 

Here, E^^iso is the isotropic gamma ray energy, di is the luminosity distance as derived with 
the concordance cosmology {fl\ = 0.7, Q^n = 0.3, Hq = 74 km/s/Mpc), Sboio is the bolometric 
fluence (the fiuence over the burst rest frame 1-10,000 keV range), and z is the redshift of 
the burst. The quantity has been called the 'energy ratio' for the Amati relation (e.g. 

•^bolo 

Band & Preece 2005). The left side of equation 3 uses only directly observable quantities 
(albeit, they are model dependent), while the right side is only a function of distance. As 
the distance rises, dj^ gets larger and (1 + z)~^'^^ gets smaller, which gives a maximum value 
for the right side. When the concordance cosmology is used, the function peaks at 2; ~ 3.6. 
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Specifically, the right side of the equation cannot exceed 1.13 x lO^keV^'^'* erg ^ cm^ and, 
therefore, 

JT2.04 

-P^t^ < 1.13 X lO^keV^ erg-^ cm^. (4) 

^bolo 

This becomes a simple way to test the Amati relation even for bursts without redshifts. 
Similarly, for the Ghirlanda relation, 

= (1.35 X 10^' erg keV-^'^^') [Ep (1 + z)]''^' (5) 

7-1 '^'^(^T.^bolO'^beam /«\ 

= TT-z 

Smo ^ (1 + z)2-43(l.35 X 1047 erg keV-i-^s) ^ ^ 

The beaming factor, Ebeam, is defined as (1 — cos6'jet), where 9 jet is the opening angle of the 
jet of the burst. The right hand side has a maximum value at Zmax = 12.6 with a value of 
2.7 X 10^° keV^'^^ erg~^ cm^ for Fheam = 1- Thus, the Ghirlanda relation forces the limit. 



r;1.43 

-ipl^ < 2.7 X 10^° keV^-^^ erg-i cm^ (8) 



bolo 

So we have a simple observational test for compliance with the Ghirlanda relation. We also 
have reproduced the result that the 'energy ratio' for the Ghirlanda differs from the Amati 
relation (e.g. Band & Preece 2005). 

At first glance, it might appear that this it is being overly generous to apply a beaming 
factor of Fbeam = 1 to calculate the limit of the Ghirlanda relation. The whole point of 
applying such a beaming factor is to give an illustration of the lowest value the limit can 
have. If we were apply some sort of average beaming factor, the value of the Ghirlanda limit 
would increase, resulting in many more rejected bursts. As an example, let us use a typical 
jet angle, 9jet = 8.5°, which corresponds to a beaming factor, Fbeam = 0.01. In this case, the 
Ghirlanda limit would increase by a power of two to ~ 2.7 x lO^keV^'^^ erg~^ cm^. There 
are two reasons, however, for choosing to keep the Ghirlanda limit with a -Fbeam = 1-0. The 
first is the simple mathematical statement that we are looking for the maximum value for 

^1.43 

The second, and more physical explanation, is that the Ghirlanda relation accounts 
for all beaming factors, and thus the limit should given as when Fbeam = 1-0. Indeed, as 
we will see, the Amati relation can be thought of as a sort of Ghirlanda relation that was 
constructed with some average beaming factor of bursts. 

Nakar & Piran (2005) analyzed 751 BATSE bursts, finding 48% of the bursts to violate 
the Amati limit (equation H]). They used this result to declare the Amati relation as not 



-5 - 



reliable. Schaefer &: CoUazzi (2007) later showed that this fraction was what is expected 
for this relation due to simple scattering effects. In the ideal case with no measurement 
uncertainties, there are no violators, but as soon as noise is introduced, some of the points 
(particularly those near the limit) would become violators. Of course, there would be an 
equal number of bursts that would go up (and thus above the limit) to those who go down 
(and below the limit). With the function on the right side of equation 3 being near its max- 
imum value (when compared to the uncertainties in Ep^ak and the scatter in the luminosity 
relations) for the redshift of most bursts (1 < z < 6), this means that about half of the GRBs 
will be apparent violators even if the Amati relation is correct. So the real finding is that 
the differences between the limit and the equation are very small given the known scatter 
(Schaefer & CoUazzi 2007). Thus, the original test by Nakar and Piran actually confirmed 
the Amati relation. 

Schaefer & CoUazzi (2007) extended this test to 69 bursts with known redshifts from 
many satellites. The result was that the Amati relation had an expected amount of violators, 
44%. The paper goes on to show that this test could be generalized for use to test all of 
the luminosity relations. In most cases, this resulted in no maximum, or at the very least 
no maximum out to very high redshift (~ 20). All luminosity relations were shown to have 
either no bursts failing the test (like the lag- luminosity relation), or a number of failures that 
is acceptable given measurement errors) . The end conclusion is that all accepted luminosity 
relations of the time had passed the Nakar and Piran test. 

A complication arises in Band & Prcccc (2005), which examined (largely) the same 
BATSE bursts as Nakar and Piran without known redshifts, finding a >80% failure rate 
for the Amati relation. This result was later confirmed by Goldstein et al. (2010). In 
particular, the spread of points was greatly over the limit, and this cannot be attributable to 
measurement error or scatter in the luminosity relation. With the differences in the violator 
fraction and the interpretation, we have a core dilemma for the Amati relation, and the 
understanding of these differences is the core of this paper. 

Another criticism of the Amati relation is that the Amati relation is both dependent on 
the sateUite and that it arises from selection effects (Butler et al. 2007). Butler et al. (2007) 
pointed out an apparent shift in the Amati relation between Swift and pre-Swift data sets. 
Their claim that selection effects will produce the Amati relation were never substantiated 
by any analysis, examples, or derivations, and the cause of the selection effects was never 
identified. These claims have scared off some workers from using any luminosity relations 
(e.g. Bromm & Loeb 2007); however, it has been widely rejected for a variety of strong 
causes (e.g. Cabrera et al. 2007; Amati et al. 2009; Ghirlanda et al. 2009; Krimm et al. 
2009; Nava et al. 2009; Xiao & Schaefer 2009) and most recently the authors themselves 
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recanted their previous findings (Butler et al. 2009). Nevertheless, we are here reconsidering 
the cause and effects for the basic claims of satellite-to-satellite differences. 

In this paper, we first start by presenting and explaining the Nakar and Piran test, 
which following Band & Preece (2005) we extend by considering bursts in a plot of their 
Sboio versus Epeak,obs- In addition, we explain why a certain amount of violators are expected, 
and what the observed distributions of bursts tell us to expect. We follow this by showing 
gathered data from various detectors and providing a comprehensive examination of how 
each detector's data performs under the Nakar and Piran test. Following this, we provide 
an explanation for why the vast majority of the data sets have too many violators of the 
Amati limit, and therefore the Amati relation is not good as a luminosity relation. Finally, 
we examine several sources of systematic offsets that are actually the cause of the Amati 
relation in the first place, which only further condemns the Amati relation's usefulness. 



2. The Shoio — Epeak Diagram and the Amati Relation 

One way to visualize bursts under the Nakar and Piran test is to plot them on a Sboio — 
Epeak,obs diagram. In doing this, we can not only easily see how a certain group of bursts 
fares on the Nakar and Piran test, but we also determine if there is a systematic offset 
between different detectors. As an example, we can determine whether different detectors 
are pre-disposed towards different regions on the diagram. In addition to plotting points 
for individual bursts, we can also plot the Amati limit (equation H]) and the Ghirlanda limit 
(equation |8]) for easier visualization of where the limits lie. Figure [T] shows the basic idea 
behind the plots, with three zones for whether the burst violates no limit, the Amati limit 
only, or both limits. 

To illustrate the Amati limit. Figure [2] presents our Monte Carlo simulation of 1000 
bursts where the Amati relation is adopted. There are no measurement errors, no selection 
effects for satellite detectors, and the burst luminosity and distance distributions are a rea- 
sonable model of the real Universe. Our simulation for each burst starts with the random 
selection of Epeak,obs as based on a log-normal distribution like in Mallozzi et al. (1995). 
In addition, the redshift of the burst is randomly selected from a log normal distribution 
centered at z=2 with a standard deviation of 1. While this is not exactly the same as the 
distributions seen in cosmological models, it is similar enough, and the exact shape does 
not effect the conclusion. With these two values, we first find the intrinsic Epeak by simply 
applying the redshift correction: Epeak = -Epeafc,o6s(l + z). We then use the Amati relation 
to derive E^^iso, and use equation 2 to get the observed Sboio- As such, the figure shows 
a realistic distribution, or at least for no measurement uncertainties. In the figure, we see 
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that there are no violators (i.e., bursts appearing below the Amati limit), with most bursts 
appearing close to the limit line. This figure is a central illustration of the Nakar & Piran 
test, which we will extend in this paper. 

If we allow for ordinary scatter caused by measurement errors in Epeak,obs and Shoio, then 
the tight scatter in Figure |2] is lost. This is shown in Figure [3l where suddenly somewhat 
less than half of the bursts become violators. For this simulation, we assumed that the 
measurement errors have a log-normal distribution with a one-sigma width of 0.25 (CoUazzi 
et al. 2011). The exact fraction of violators will depend on the size of the observational 
scatter. In this realistic simulation, ~40% of the bursts are below the Amati limit line. The 
point of this figure is that normal and expected observational measurement errors will lead 
to nearly half the bursts being apparent violators. Importantly, this scatter does not explain 
the high violator rates reported by Band & Preece (2005) and Goldstein et al. (2010). This 
discrepancy is the main topic of this paper. 

For comparison, we can also consider how the Sboio — Epeak,obs diagram would look if 
neither the Amati or Ghirlanda relations were valid. For this, we have constructed another 
Monte Carlo simulation (see Figure H]). As in Figure 2, we have assumed no measurement 
errors, no selection by satellite detectors, and we have adopted realistic luminosity and 
distance distributions, but we have made no constraints from either the Amati or Ghirlanda 
relations. We start by selecting burst distances and energies in the 100-500 keV such that 
they reproduce the observed log(A^) — log(P) curves for BATSE (Fenimore et al. 1993, 
Fishman & Meegan 1995). We then generate -Epeafc based on a log normal distribution with 
some loose connection to the brightness of the burst (as seen in Mallozzi et al. 1995). We 
then apply a bolometric correction with (a = —1.0 and /3 = —2.0). The result is an in 
illustration of the intrinsic distribution of bursts on the sky. Our simulation of 10,000 bursts 
has approximate edges at 20 and 3000 keV, plus lower and upper edges simply where we cut 
off the log(A^) — log(P) curve. The key point is that Figures 2 and 4 are greatly different, 
because low-fluence bursts will dominate unless some law/correlation forces these low-fluence 
events to have \ow-Epeak,obs- So we have two extreme cases that produce greatly different 
distributions in the S^oio — Epeak,obs diagram. 

Both Figures 2 and 4 are for the intrinsic distributions of GRBs in a realistic case with no 
effects of detector thresholds or measurement uncertainties. From a comparison of Figures 
2 and 3, we see that realistic measurement errors will substantially smear the underlying 
distribution. Detector thresholds will also force a fuzzy cutoff roughly running along some 
horizontal curve. For a detector with a high threshold, the many violators in Figure 4 will 
never be detected and the violator fraction might appear acceptable. For a detector with a 
low threshold, we should be able to easily determine whether the Amati relation is valid. 
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3. Generalizing the Test to Many Detectors 

So far, the Nakar & Piran test has only been applied to BATSE bursts (Nakar & Piran 
2005, Band & Preece 2005; Goldstein et al. 2010) and to a collection of bursts with redshifts 
as detected by a range of many satellites (Amati 2002; 2006; Schaefer & CoUazzi 2007). But 
this test can be extended to many satellites, because all that is needed are values of Sboio and 
Epeak,obs, with both of thcse being commonly reported for many bursts. The essence is to find 
the fraction of violators, ^, for each sample. We will also keep track of the quantity (log -^f^) 
for each sample, as this can be directly compared to 1.13 x 10^ (in units of keV^"^ erg~^ cm^) 
so as to test the Amati limit (cf. equation 4). We will keep track of these statistics for both 
bursts with spectroscopically determined redshifts (Greiner 2010) and those with no known 
redshift. The samples and their statistics are presented in Table 1, with some discussion in 
Sections 3.1 to 3.8. The Sboio — Epeak,obs diagrams for each sample are presented in Figures 
5-13. 

3.1. Amati et al. (2006) 

We start by using the compilation of data from Tables 1 and 2 of Amati et al. (2006). 
These bursts all have redshifts and are from BeppoSAX Konus, HETE, BATSE, and Swift. 
Using the Amati relation, we calculate the Sboio from the given data. We exclude bursts 
050315, 050824, 050904, 981226, 000214 and 030723 because only limits to ^peafc,o6s are 
provided, and therefore are not useful to us. We also had to exclude burst 980329 because 
a redshift range is given, and we could therefore could not get a accurate measurement of 
{d^) for converting Eiso into Sboio. The results are in Figure O 

This sample of GRBs was largely the same as used by Amati (2002) to discover and 
calibrate the Amati relation, so it is no surprise that the bursts are spread out along the 
Amati limit line. The violator fraction is ^ = 34%, which is as expected given the usual 
scatter due to measurement errors. The sample has (log = 8.90 which is close to the 

^bolo 

£;2.04 

limit of log(1.13 x 10^) = 9.05 (in appropriate units). The RMS scatter of log is 0.56, 
which is a measure of how tight the Amati relation is for the sample. 

3.2. Schaefer (2007) 

Next, we visit another compilation set, this time from Schaefer (2007). This data set 
also takes its burst sample from a variety of different detectors: Konus, BATSE, Beppo-SAX, 
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HETE, and Swift bursts were included for this sample. While the paper studies 69 GRBs 
with known redshift, only 27 have the bolo metric fiuence reported and thus are the ones that 
we use (see Figure ED- 

Of the 27 bursts, 11 fail the Nakar and Piran test (^=41%). This is an expected failure 
rate for the Amati relation and is in agreement with previous analysis on this data (see 
Schaefer and Collazzi 2007). The average value for the log of (log '"="'' ) is 8.95±0.57. So 
this data is similar to the previous data set (which is not surprising as they share some of 
the same bursts). 

The two samples which contain exclusively bursts with known redshifts both agree well 
with the Amati relation, and this has long been the primary justification for accepting 
the Amati relation as a physical relation for GRBs. However, other samples (see below) 
do not agree with the Amati limit, and this suggests that bursts with redshifts might be 
a significantly different sample from those without redshift. This is the reason why we 
distinguish bursts with and without redshifts in Table 1 and in the Figures. 

3.3. BATSE data 

Our data for BATSE is a part of the upcoming 5B catalog (Goldstein et al. 2011). We 
present the values of Epeak and Sboio for the most statistically preferred fitting model, GPL 
or Band. We also use bursts for which we have 40% relative errors or better. After applying 
these selection criteria, we are left with 1654 bursts, which are plotted in Figure [71 In the 
figure, we also introduce two new curved lines which represent illustrative thresholds for the 
trigger (dotted) and the ability to detect Ep^ak (dot-dashed). These are explained in detail 
in section 4. 

We find that the BATSE bursts fail at an extreme rate, with 93% violators. In addition, 
the BATSE bursts cover a large region of the disallowed zone, with very few bursts above the 

^2.04 

limit. We find that (log J"""'' ) to have a value of 10.18±0.88. The failure rate is consistent 
with that observed in the past of BATSE bursts in previous works. The spread of BATSE 
bursts is so large, it hints that most any future changes to the Amati relation (e.g. as more 
Swift bursts are detected with redshifts) will result in a high failure rate. 
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3.4. 



HETE data 



Our sample of HETE bursts comes from Sakamoto et al. (2005). The quoted values for 
Sboio were covering the 2-400 keV range, and had to be converted into bolometric fluences. 
This was done by using the given parameters for the spectral model (Band or Cut-off Power 
Law (CPL)) to extrapolate a bolometric correction. This was generally a small correction, 
while even the large corrections are still small compared to scatter in Figure El The quoted 
error bars for both Epeak and Sboio are given for the 90% level, so they had to be converted 
into standard one sigma confidence level error bars. The figure also has illustrative lines 
to represent trigger thresholds (dotted line) and the Ep^ak detection threshold (dot-dashed 
line). Again, these will be explained in detail in section 4. 



The HETE bursts have ^ = 33% and (log -^^) similar to that of the original Amati 



sample, so it appears that these bursts are consistent with the Amati relation. Nevertheless, 
the scatter apparent in Figure |H] is so large that there is little utility in applying the Amati 

^2.04 

relation to these bursts. The difference is insignificant for (log J"'"'' ) between bursts with 



and without redshifts. Of all the single-satellite data sets that we consider, the HETE data 
is the only one that apparently obeys the Amati relation, although its large scatter limits 
its usefulness. 



For the Swift data, we use the catalog in Butler et al. (2007). We use their Epeak,obs as 
derived from frequentist statistics as it is the most common approach to finding Ep^ak- (Their 
Bayesian values were made with unreasonable priors that significantly skew the results.) The 
Swift burst detector only goes up to 150 keV, so the reported values of Epeak,obs are almost all 
lower than 200 keV. We have adopted their bolometric fluences and we have converted their 
non-standard 90% error bars into standard one-sigma error bars. Figure [9] plots the results; 
bursts without known redshifts are represented as empty circles and bursts with known 
redshift are represented by a filled diamond. This will be true for all future plots. This is the 
last of the three plots in which we have plotted illustrative lines representing trigger (dashed 
line) and Epeak detection thresholds (dot-dashed line), which again, are detailed in section 
4. 

The Swift bursts violate the Amati limit at a rate of 76% to 82%. That is, the Amati 
relation does not work for Swift. This result is not the result of small number statistics, and 
we can see from the distribution that the disagreement is highly significant. This is another 
version of the same conclusion first reported by Butler et al. (2007). 




3.5. Swift 
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When first confronted with the discrepancy that bursts with redshifts agreed with the 
Amati relation (see Figures |5] and |6]) while bursts without redshifts disagreed with the Amati 
relation (Band & Preece 2005), our initial thought was that the bursts with redshifts might be 
somehow selected from a separate population for which the Amati relation applied. However, 
with this large sample of well-measured bursts from Swift, the distributions of bursts with 
and without redshifts is essentially identical. Thus, we have a proof that the success or 
failure of the Amati relation does not depend on some selection effect that correlates with 
the measuring of spectroscopic redshifts. Another thing to remember is that since Swift 
bursts are the bursts that account for a majority of the bursts with known redshifts, there is 
a built-in selection effect that will eventually develop that will bias future iterations of the 
Amati limit towards the area Swift bursts cover. 

3.6. Suzaku 

Suzaku data for long GRBs are available through the GCN circulars (|http: / / gcn.gsfc.nasa.gov ) . 
Typically, the reported fluence covers the 100 keV to the 1 MeV range, so we apply a bolo- 
metric correction based on the reported spectral fit. A typical bolometric correction value is 
a factor of ~1.7. The -Epeafc,o6s and fiuence values reported in the circulars are preliminary 
and made soon after the burst, yet any likely changes to get to the final best fits are greatly 
smaller than the scatter shown in Figure [10] and are thus not important. The Suzaku bursts 
all have Epeak,obs > 200 keV, a result of the relatively high energy range of sensitivity of the 
detector. 

Most of the Suzaku bursts violate the Amati limit, usually by a large factor (^ ~, 94%) 
and the small fraction that are not violators are very close to the limit (Figure [TU]) . Therefore, 
the Amati relation does not work for Suzaku bursts. This is true for both bursts with redshifts 
and without. 



3.7. Swift-Suzaku 

Krimm et al. (2009) presented a catalog of bursts for which there was both Swift and 
Suzaku data. The expanded energy range gave better fits to the spectra, with Swift covering 
the lower energies and Suzaku covering the higher energies. With the joint spectral fits over 
a very wide range of photon energies, the sample has a wide range of Epeak,obs values from 30 
keV to 2000 keV. The catalog lists the best fit for the three major spectral models (power 
law, power law with an exponential cutoff, and Band model) for the majority of the listed 
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bursts, while we use the Epeak,obs values from just the Band function. 

The joint sample largely stretches between the Amati limit line and the Ghirlanda limit 
line (Figure [TT]) . The fraction of violators is 86% for the 28 bursts with spectroscopic redshifts 
and 74% for the 38 bursts without redshifts. Again, the Amati relation fails, and there is 
no significant difference related to whether the burst has a spectroscopic redshift or not. No 
burst significantly violates the Ghirlanda limit. 

3.8. Konus 

The GRB detectors on the Wind satellite (Aptekar et al., 1995) are long-running instru- 
ments with a stable background that has measured many bursts, with the fluence and Epeak,obs 
values promptly reported in the GCN circulars. These reported value are preliminary, with 
no final analysis having been published, but any plausible errors due to the preliminary 
nature of the report are greatly smaller than the observed scatter in the Sboio — Epeak,obs dia- 
gram (see Figure [T2|) . We find a total of 97 bursts, 33 of which have associated spectroscopic 
redshifts, with reported fiuences and Epeak,obs for the entire burst interval. We applied a 
bolometric correction factor, based on the given spectral fits. This factor was typically very 
small, due to the large range the reported Konus fiuences usually cover. 

The distribution of the Konus bursts in the Sboio — Epeak,obs diagram has a flat lower 
cutoff, likely due the trigger threshold (although this cutoff is higher the reported trigger 
threshold in Aptekar et al. (1995) of ~ 5 x 10~''ergcm~^). The distribution also shows 
a fairly high upper limit on Epeak,obs due to the sensitivity of the detectors to high photon 
energies. From 22% to 27% of the bursts are above the Amati limit line, while all the bursts 
are above the Ghirlanda limit line. The Amati relation fails for the Konus bursts. Again, 
the bursts with redshifts are distributed identically to those without, so there is no apparent 
selection effect based on spectroscopic redshifts. 

3.9. Beppo-SAX 

Guidorzi et al. (2010) provides a large catalog of both Ep^ak and S from the Beppo-SAX 
GRBM. In the catalog, data for the brightest 185 bursts are given; for which we are able 
to use 129. Of the useable 129 bursts, we have only 10 bursts with spectroscopic redshifts. 
The provided Shoio were over the 40-700 keV range. We apply the same type of bolometric 
correction as before, using the provided spectral indices for the GPL used. This correction 
is typically small, with a typical correction value of 1.5. 
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While it is impossible to make a strong statement about the shape of the distribution 
of Beppo-^KX. bursts with only the bright bursts, there is still an important result from the 
data. The data is plotted in Figure [T^ We find that even among the brightest bursts, 90% of 
bursts with redshifts and 90% of bursts without redshifts are violators. What is particularly 
provocative about this result is that these are the brightest bursts, and thus the most likely 
to not be violators. It is unlikely that there are a significant number of 'missing' bursts that 
would be non-violators. Any such burst would have to be both bright and have a low Ep^ak 
while still being sim enough to be missed in the bright burst catalog. Finally, we provide 
information as to the average energy ratio of these bursts, but we once again stress that 
these are only the brightest bursts, and therefore these values should be taken with caution. 
We therefore feel confident in saying that the Amati relation fails for Beppo-^KX. bursts, 
although this statement is not as strong as it is for other detectors because of the sample 
used. 



3.10. Overview of Results 

Previously, Butler et al. (2007) had pointed out that the normalization constant for the 
Amati relation was slightly different depending on whether Swift or pre- Swift bursts were 
used, and this is like noting that (log has changed. Previously, Band & Preece (2005) 

and Goldstein et al. (2010) pointed out that > 80% of BATSE bursts violate the Amati limit. 
In this section, we have generalized these analyses, both to looking at many GRB detector 
instruments and to looking at the two dimensional distribution in the Sboio—Epeak,obs diagram. 

All of these data sets give consistent conclusions: (1) The distribution of bursts in the 
Sboio — Epeak,obs diagram varies significantly and greatly from satelhte-to-satellite. (2) The 
only data sets to pass the generalized Nakar & Piran test for the Amati relation are the 
early heterogeneous sample of bursts with measured spectroscopic redshifts. (3) The bursts 
detected by BATSE, Swift, Suzaku, and Konus all have a high fraction > 70%) of bursts 
which violate the Amati limit, with the violations being highly significant and by large 
factors. That is, the Amati relation fails for bursts from these four satellites. (4) The Amati 
limit is satisfied for the HETE bursts, to the extent that the violator fraction is consistent the 
Amati relation plus normal observational scatter, however, the scatter in the Sboio — Epeak,obs 
diagram is so large that we conclude that the Amati relation does not satisfactorily apply 
to the HETE data. (5) We find that no bursts, from any satellite, significantly violate 
the Ghirlanda limit. (6) These conclusions are true whether we examine only bursts with 
spectroscopic redshifts or without redshifts. 

The normalization factor for the Amati relation will scale closely with (log ) . These 
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values are listed in Table 1, and can be taken as an intercept with some dispersion. That 
is, these values can possibly be seen as a means of describing a sort of Amati relation to 
each detector. However, this is not a useful interpretation as the dispersion is often as big 
as the population itself, rending such an interpretation nearly meaningless. If anything, 
the results from Table 1 demonstrate that the Amati relation must vary from detector-to- 
detector by over an order of magnitude. With the bursts seen in the sky not depending on 
the satellite, the large variations in the Amati relation from detector-to-detector imply that 
there must be some selection effect which biases the visible bursts, with these biases being 
instrument-specific. 

Every burst detector has a substantially different distribution of bursts in Figures [7] - 
[121 Since the population of bursts that appear in the skies above the Earth does not change 
with the satellite, so the large changes from detector-to-detector can only be due to some 
selection effect where bursts in various regions of the Sboio — Epeak diagram are not selected. 
The next section will investigate and identify these selection effects that create the Amati 
relation. 



4. The Amati Relation Comes from a Combination of Selection Effects 

The distributions of bursts in the Sboio — Epeak,obs diagram are caused by a variety of 
effects. Some of these effects are caused by detector limitations that prevent a burst from 
appearing in some parts of the Sboio — Epeak,obs diagram (Sections 4.1 and 4.2), while other 
effects make for rare bursts in other regions of the Sboio — Epeak,obs diagram (Sections 4.3 and 
4.4). The combination of these effects will produce the observed distributions (Section 4.5). 
For some detectors, the selection effects will force the observed bursts to follow a roughly 
diagonal region (with wide scatter) that will appear as the Amati relation (Section 4.6). 

4.1. Trigger Thresholds 

The best known selection effect is the detector trigger threshold. For example, a burst 
would trigger BATSE only if it was produced a peak flux (in a 0.064, 0.256, or 1.024 second 
time bin) brighter than 5.5-a above background in at least two detectors over the 50-300 keV 
energy range. Other satellites have more complex trigger algorithms (for example, GBM has 
overlapping triggers), but they all come down to the same essentials. The trigger threshold 
depends on the Epeak,obs, the spectral energy range of the trigger, the background flux, and 
the effective area of the detector. The triggers operate off the peak flux (Pmax), so the 
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limiting fluence will depend on the effective duration (Sboio/Pmax), which can vary widely 
from burst to burst. Thus, the limit due to trigger thresholds will be 'fuzzy', with no sharp 
edge but rather a gradient as Sboio is reduced. Approximately, the trigger threshold will 
produce a horizontal cutoff at the bottom of the Sboio — Epeak,obs diagram. 

In principle, the exact trigger thresholds can be calculated for every detector and burst. 
In practice, the conditions {Epeak,obs, background flux, incidence angle, burst light curve) vary 
greatly from burst to burst, creating substantial scatter in the thresholds. For this paper, 
we do not need an accurate distribution for the Sboio threshold, so instead we calculate the 
typical Sboio threshold as a function of Epeak,obs for average conditions. In particular, we 
adopt an average spectral shape as the Band function (Band et al. 1993) with a low-energy 
power law index of -1.0 and a high-energy power law index of -2.0. We also take the effective 
duration of the peak {Sboio/ Pmax) such that it fits the observed distribution of the detectors. 
For each detector, we take its trigger energy range, face-on effective area, and the average 
background flux. The formalism and many of the input parameters were taken from Band 
(2003). The result is a lower limit in the Sboio — Epeak,obs diagram, as displayed in Figure [7] 
for BATSE, Figure [S] for HETE, and Figure [H] for Swift. We do not have enough information 
to calculate trigger thresholds for some satellites, but the threshold is usually fairly obvious 
(e.g.. Figure [12] has a nearly flat and moderately sharp lower limit to Sboio)- These thresholds 
are not sharp, so bursts can easily appear somewhat below the threshold. Indeed, by varying 
the input conditions somewhat, the trigger threshold lines can be translated up and down 
substantially. 



4.2. Threshold For Measuring Epeak,obs 

A second detector selection effect is that the burst must have enough photons recorded 
for the analyst to be able to determine the Epeak,obs value. This will depend on both Sboio 
and Epeaks,obs as well as the detector properties. For example, a burst just above the trigger 
threshold will have just enough photons to be detected but not enough photons to allow 
any constraints on the Epeak,obs value, so this burst will not be included in a sample for 
plotting on the Sboio — Epeak,obs diagram. For another example, consider a burst with Epeak,obs 
at the upper edge of the measured spectral range for a detector, such that a very bright 
burst will have a well-measured turnover that accurately deflnes the fltted Epeak,obs value, 
whereas a fainter burst will have poor photon statistics near the turnover in the spectrum 
and the Epeak,obs value will remain unmeasured and the burst will not be included in any of 
our samples. 

In general, for a given Epeak,obs, there will be some lower limit on Sboio, below which 
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there will be too few photons to measure Epeak,obs- As Epeak,obs moves to higher energies, 
the limit on Sboio will sharply increase. The result will roughly be a diagonal line across the 
Sboio — Epeak,obs diagram, from lower left to upper right, with any burst below that line not 
having a measured -Epeafc.ofes and not appearing in any sample of bursts in Section 3. 

We have made calculations of this threshold curve for BATSE, HETE, and Swift. To 
do this, in a Monte Carlo sense, we constructed many simulated bursts over each detector's 
spectral range for many values of Epeak,obs where the normalization and error bars of the 
spectra were determined by the burst fluence. These spectra were also for each spectrum, 
we then fitted both a power law times an exponential model (with the calculated Epeak,obs 
value) and a simple power law model. If the chi-square values for the two fits differed by 
more than 15.0 (so that the model with the peak was a sufficiently good improvement on 
power law model given the extra degree of freedom), then we took the Sboio for the burst to 
be above the threshold. By varying the Epeak,obs, we were able to determine the threshold 
for measurement as a curve in the Sboio — Epeak,obs diagram. As these fines are merely for 
illustration, we do not use the DRM for these simulations. For BATSE, HETE, and Swift, 
our calculated thresholds are presented as curves in Figures [3, [HI and M 



4.3. The Epeak,obs Distribution 

Amongst bursts appearing in the skies, the Epeak,obs distribution is not flat, but rather 
bursts appear with a roughly log-normal distribution of Epeak,obs- For bright bursts, the 
mean value is 335 keV, with the FWHM stretching from roughly 150-700 keV (Mallozzi et 
al. 1995). This mean value shifts signiflcantly as the bursts get dimmer, being 175 keV just 
above the BATSE trigger threshold (Mallozzi et al. 1995). The so-called 'X-ray Flashes' are 
simply bursts in the low-energy tail of the distribution (Kippen et al. 2004; Sakamoto et al. 
2005; Pelangeon 2008). The existence of this single peak in the Epeak,obs histogram is highly 
signiflcant and not from any instrumental or selection effect (Brainerd et al. 1999). In all, 
most bursts are between 100-700 keV, and bursts <30 keV or >1000 keV are rare. This will 
directly translate to unpopulated regions of the Sboio — Epeak,obs diagram. A direct simulation 
of this distribution is given in Figure HI 

The Epeak,obs distribution wfil cause deflnite but gradiated cutoffs in the Sboio — Epeak,obs 
diagram. These cutoffs will be nearly vertical. The drop in the average Ep^ak,obs will make 
the cutoff on the right have a slope down to the lower left. 
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4.4. The Sboio Distribution 

Unsurprisingly, bright bursts are rare, while faint bursts are more frequent. The dis- 
tribution of burst fiuences is traditionally represented by the log A^(> P) — logP curve, for 
which the best observations come from the BATSE catalog (Fishman et al. 1994; Paciesas 
et al. 1999). For bright bursts, the slop of the curve is nearly the ideal -|. The slope flattens 
out for faint bursts, approaching -0.7. In the Shoio — Epeak,obs diagram, the density of bursts 
falls off drastically from bottom to top (see Figure H]). 

4.5. The Effects in Combination 

The intrinsic distribution of bursts in the 6*60/0 ~ Epeak,obs diagram is determined by 
the Epeak,obs log-uormal distribution that changes with Sboio (Mallozzi et al. 1995) and by 
the logA^(> P) — logP distribution (Fishman et al. 1994). With these two effects, the 
burst density across the diagram is displayed in Figure HI Together, the two effects produce 
contour lines of burst density in the diagram, and two such curves are displayed in Figure [HI 
The combined effects make for the upper-left corner of the Sboio — Epeak,obs diagram, simply 
because bursts up there are doubly rare (both low in Epeak,obs and very bright). This means 
that there are few bursts significantly above (and to the left) of the Amati limit line. That 
is, there is a natural cutoff on the top side of the Amati relation. 

The detector selection effects then operate on the natural distribution. The well known 
trigger threshold is actually below the threshold for measuring an Epeak,obs value, so only the 
later selection effect is really operating. This selection effect cuts on a sort of diagonal from 
lower-left to upper- right, and its position depends greatly on the detector sensitivity and 
energy band for the trigger. For a relatively poor sensitivity and a trigger energy band that 
effectively does not get much above a few hundred keV, the threshold will be quite high. 
Indeed, for many detectors, the threshold will be just below the Amati limit line (Figure 
[T^ . so there will be few bursts significantly below the Amati limit line. That is, there is a 
selection effect from the intrinsic distribution of bursts such that there is a natural cutoff 
below the Amati limit. 

For some detectors, we can see that the Amati relation is a natural and expected con- 
sequence of the intrinsic burst distribution combined with normal detector selection effects. 
This is illustrated in Figure [TH where the allowed region is confined to an area along the Am- 
ati limit line. From Figure El we know that bursts along the Amati limit in the Sboio — E peak, obs 
diagram will then imply a relation close to the Amati relation. Thus, the natural distribution 
of bursts makes for bursts above the Amati limit (i.e., the very-bright \ow-Epeak,obs bursts) 
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to be rare, while the detector selection effects makes for bursts below the Amati limit (i.e., 
the faint high- Epeak,obs to be too faint to have a measured Epeak,obs- With the only bursts 
remaining being close to the Amati limit line, a relation like the Amati relation would be 
apparent. Thus, we conclude that the Amati relation is simply a result of selection effects 
and there is no physical basis. 

For the original bursts used to define the Amati relation (Amati 2003), an additional 
selection effect is operating, in that the burst must also have a measured spectroscopic 
redshift for inclusion in the sample used for calibration. The selection effects for measuring 
a redshift are complex. There is certainly a selection based on redshift, with the cause being 
that more distant bursts are fainter, hence less likely to have a visible optical transient or 
host galaxy bright enough to get lines in a spectrum. An additional effect on redshift relates 
to the availability of spectral lines in the optical band. The efficiency of measuring redshift 
as a function of redshift has been quantified in Xiao & Schaefer (2011), with this effect 
being roughly an order of magnitude between z = b and nearby bursts. The efficiency of 
measuring redshifts also presumably depends on Sboio which will roughly scale as the burst 
brightness in the optical band. The redshift measurement efficiency also is time dependent, 
as optical follow-up strategies and capabilities change within our community. Thus, Swift 
bursts started out with an average redshift of 2.8 in the first year after its launch (Jakobsson et 
al. 2006), while the average redshift has steadily declined to 2.1 over the last year (Jakobsson 
et al. 2008). The reason for this shift is unknown, but it must come from overall follow-up 
practices in our community. The bursts with redshift used in the original calibration of the 
Amati relation have an average redshift of 1.5, indicating that the effective threshold for this 
sample is quite high. 

The distribution of bursts in the S^oio — Epeak,obs diagram will depend greatly on the 
detector. The threshold for measuring Epeak,obs varies substantially detector to detector. For 
example BATSE has a low threshold while Konus has a high threshold. The shape of the 
threshold (as a function of Epeak,obs) also varies, from a flat bottom for Konus due to its 
sensitivity to high energy, to the up-sloping threshold for Swift due to its lack of high energy 
sensitivity and even more exaggerated in HETE with it's small area. The ability to measure 
Epeak,obs dcpcuds Critically on the energy range of the spectra. The Konus detectors have 
a very wide range of spectral energy resulting in a wide range of measured Epeak,obs values, 
the Swift detectors cutoff around a few hundred keV, while the Suzaku detectors can only 
record Epeak,obs ^ 200 keV. The combination of these selection effects makes the distribution 
of bursts different for each detector, and accounts for the wide range of distributions seen in 
Figures 

Still, the issue has been raised in recent tests (e.g. Ghirlanda 2011) that what we are 
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seeing in these failures is just the scatter about a relation which is ever changing with new 
bursts every day. The primary argument is that the Nakar and Piran test limit should be 
formulated from the 3-sigma line about the model, instead of the model line itself. As a 
result, the Amati limit would be considerably higher. There are a variety of problems with 
this argument. The first of which being that there is already an allowance made for the 
Amati relation to have up to 40% violators and not be considered as failing for the data 
set. Therefore, the scatters are already being accounted for, and it is overkill to use such 
a generous limit to perform the test. If the test is done in this manner, no longer can 
allowances be made for any violators (or, more precisely, there needs to be less than 0.3% 
violators). Even by the groups own tests, there are violators on the order of a few percent, 
depending on the test. This is an unacceptable violator rate considering they are violating 
a limit from the three-sigma deviation from the model. Finally, another question that arises 
is that the bursts wc sec all seem to be biased in one direction. If we were seeing the result 
of measurement scatter about the Amati relation, we should expect to sec an equal fraction 
of bursts well above the limit line. Instead, we see that for almost all data sets, the bursts 
are systematically in one direction from the limit. 

The Amati relation will certainly see improvement in these tests in the future. With 
increasing number of Swift bursts with spectroscopic redshifts, it will undoubtedly eventually 
lie right in the middle of the Swift data set. Even then, the Amati relation will be failing 
for our best data sample, the BATSE data. Our argument is that there are undeniable 
systematic effects at play that are causing the Amati relation, and therefore even these 
'improvements' would be fairly meaningless as we would still see systematic differences in 
where bursts are observed in the diagrams. Therefore, the Amati relation is simply not good 
for making any kind of predictions, cosmological or otherwise. 

Perhaps the simplest disproof of the Amati relation is simply that the violator fraction is 
greatly too high in most data sets. And perhaps the simplest proof that the Amati relation 
is caused by selection effects is the large differences between the various Sboio — Epeak,obs 
diagrams for the many detectors. 

5. Conclusions 

The Shoio — Epeak,obs diagram has two limit lines, where bursts cannot be below that line 
if the Amati or Ghirlanda relation holds. Actually, with the fairly large total uncertainties, 

substantially larger than the simple measurement errors quoted in the literature, we can 
expect nearly half of the bursts to be scattered below the Amati limit line. So a simple 
test of the Amati relation is whether the average burst falls below the Amati limit. (This is 
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similar to the original test proposed by Nakar &; Piran, except that agreement with the Amati 
relation corresponds to about 40% violators.) We apply this test to many burst samples. 
The samples of early bursts with spectroscopic redshifts (as originally used to calibrate the 
Amati relation) pass our test, as does the sample of HETE bursts (even though the scatter 
about the Amati relation is unusably large). All other satellites have a large fraction of 
violators far below the Amati limit line. This is true whether we look at bursts with or 
without measured spectroscopic redshifts. This constitutes a proof that the Amati relation 
could possibly apply, at best, to only a small and unrecognizable fraction of GRBs. Indeed, 
the wide variations in distribution from detector to detector constitute a proof that selection 
effects must dominate the Amati relation. 

We find four selection effects restrict the distribution on all sides. The best known 
detector selection effect is the trigger threshold, which produces a roughly horizontal and 
fuzzy cutoff. A more subtle and more restrictive selection effect is that for an Epeak,obs value 
to be reported, the burst must be brighter than some threshold, with this threshold rising 
fast with increasing Epeak,obs- These two detector selection effects will cut out bursts that 
are some combination of faint and hard, with these effects changing greatly from detector to 
detector. The third and fourth selection effects operate to restrict the burst population as 
it appears in the sky. The third selection effect is that bursts have a log-normal distribution 
of Epeak,obs with the mean value shifting to lower values for faint bursts. This effect will also 
reduce the number of detectable bursts that are faint and hard. The fourth selection effect 
is that bright bursts are much rarer than faint bursts, as quantified by the usual power-law 
log[A'"(> P)] — log[P] curve. The combination of the third and fourth effects means that 
the bright and soft bursts are doubly-rare, so that the upper-left side of the Sj^oio — -Epeafe,o6s 
diagram will be empty. 

For a detector with a range of spectral sensitivity and a low detection threshold, the 
distribution in the Siyoio — Epeak,obs diagram will extend relatively low, with a large fraction 
of violators below the Amati limit (like for BATSE). For a detector with a low energy range 
of sensitivity and a low detection threshold, the cutoff will be a diagonal line just below the 
Amati limit. When combined with the paucity of bright-soft bursts in the GRB population 
(i.e., those above the Amati limit line), we have a combined selection effect that picks out 
bursts near the Amati limit. Such a burst sample would then appear to follow the Amati 
relation. Thus, the very strong selection effects for the early bursts with spectroscopic 
redshifts will create the Amati relation without any need for a physical connection between 
the Epeak,obs ^-nd Sboio- That is, the Amati relation is not real, but its appearance in some 
data sets is simply a result of various selection effects by the detectors and within the GRB 
population. 
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With these strong results, the Amati relation should clearly not be used for purposes of 
cosmology, as has been previously done by many groups. In particular, BATSE provides the 
strongest evidence for not using the Amati relation for cosmological purposes. We note that 
our groups have never used the Amati relation for any cosmological purpose (e.g., Schaefer 
2007; Xiao & Schaefer 2011). 

Our community has expected a connection between Ep^ak and energy for some time 
(e.g. Lloyd et al. 2000). With the Amati relation clearly unsuitable, the other possible 
interpretation is non-physical correlation between Ep^ak and Eiso- As demonstrated above, 
this relation is cxtrcmclv weak (e.g., sec Table 1). That is, for any given detector, some 
(log -^^\ of observed bursts can be used as a sort of intercept for such a correlation with 

^bolo 

the given dispersion. For some detectors, this empirical correlation can serve as a useful 
description of the data, and indeed this is what gave rise to the proposed Amati relation. 
However, for many detectors, this interpretation is not useful, however, because the given 
dispersions are as wide as the data itself, so any such statements regarding a correlation 
between the two would be largely meaningless. 

We emphasize that the failure of the Amati relation in no way carries any implications 
for any other GRB luminosity relation. The fault of the Amati relation can be viewed as if it 
is merely a version of the Ghirlanda relation except that the beaming correction is unknown, 
so isotropic emission was assumed. The result, however, is that the Amati relation is biasing 
itself towards some average of whatever the beaming factors of the calibrating bursts are. 
All the other GRB luminosity relations do not involve beaming corrections, and the known 
physics of the beaming is already accounted for in the physics derivations of these laws. 
The Ghirlanda relation is in essence just a conservation of energy statement, while the other 
luminosity relations (all involving the peak flux, not the fluence) just involve relativistic 
effects in the visible region of colliding jets. Indeed, most of the other GRB luminosity 
relations were predicted from the physics and then later confirmed. In all, the failure of 
the Amati relation is zero evidence for the validity of the other relations (many of which 
were confirmed predictions) and there are good physical reasons to know that they are valid 
physical laws for GRBs. 
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Table 1. Demographics of the Data Samples 



Data Set 


#w/Z 




p2.04 


H^w/o Z 


iw/o z" 


g2.04 


Ideal, no scatter 




0% 


8.95±0.29 




0% 


8.95±0.29 


Ideal, with scatter 




~40% 


9.33±0.61 




~40% 


9.33±0.61 


Amati et al. 2006 


50 


34% 


8.90±0.56 









Schacfcr 2007 


27 


41% 


8.95±0.57 









BATSE 









1654 


93% 


10.18±0.88 


HETE 


12 


33% 


8.67±0.62 


24 


54% 


9.05±0.84 


Swift 


25 


76% 


9.42±0.47 


46 


85% 


9.46±0.45 


Suzaku 


7 


100% 


9.77±0.74 


25 


92% 


10.28±0.87 


SwiftSuzaku 


28 


86% 


lO.Olibl.Ol 


38 


74% 


9.63±0.85 


Konus 


33 


73% 


9.42±0.58 


64 


78% 


9.68±0.87 


Beppo-SAX 


10 


90% 


9.36±0.39 


119 


90% 


9.51±0.39 



is the fraction of bursts that violate the 'Amati limit' of 1.13 x 10^ keV^ erg ^ cm^ for 
the Amati energy ratio E'^^^^j Shoio- 

^For the Amati relation, the quantity E'^'eakl ^bolo should never exceed 1.1 x 
10^ keV^''^'^ erg^"*^ cm^, so even with normal observational scatter a sample of GRBs should 
not have (log-^^^) > 9.04 in appropriate units. With a reasonable distribution of burst 
distances, the limit will be even smaller. So this column provides a measure of the disagree- 

g2.04 

ment with the Amati limit for a sample of bursts. The RMS scatter for log is given 
after the average value, so this can give a measure of the scatter of the distribution. 
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Epeak,Obs (keV) 



Fig. 1. — The basics of the Nakar and Piran test in graphical form. Any burst (even without 
a known redshift) can be plotted on this diagram. If the Amati relation is correct, then 
any burst must lie above the solid line (from equation Hj), although normal scatter from 
measurement error will put somewhat less than half of the bursts just below the limit line. 
If the Ghirlanda relation is correct, then any burst must lie above the dashed line (from 
equation |8]), although normal scatter from measurement error can put a small fraction of 
the bursts just below the limit line. If a burst lies below one of the lines, then it is called a 
'violator' of that relation. 
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Fig. 2. — 1,000 simulated bursts based on the Amati relation with no measurement errors. 
We assumed that the Amati relation is exact. In our Monte Carlo simulation, each burst 
had a redshift chosen randomly from a reasonable cosmological model for bursts, an Epeak,obs 
value chosen randomly from the log- normal distribution of Mallozzi et al. (1995), the burst 
energy calculated from the Amati relation, then the observed Sboio calculated from the burst 
energy and redshift. The simulated bursts are usually close to the Amati limit line, and 
there are no violators. 
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Epeak,Obs (keV) 

Fig. 3. — 1,000 simulated bursts based on the Amati relation with realistic measurement 
errors. These 1000 bursts are identical with the bursts in the previous figure, except that a 
statistical scatter has been added to the intrinsic values. For this Monte Carlo simulation, 
the measurement error is taken from a log-normal distribution where the one-sigma scatter 
in log{Epeak,obs) IS 0.25 (c.f., CoUazzi et al. 2011). With this, the fraction of burst that 
violate the Amati limit rises from zero to ~40%. 
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Fig. 4. — 10,000 simulated bursts without the Amati relation. We start by generating a 
flux in the 100-500 keV range using the BATSE log[A'"(> P)] — \og[P] relation (Fenimore et 
al. 1993, Fishman and Meegan, 1995). We then generate an Ep^ak based on a log normal 
distribution, with a dependence on the brightness of the burst (as seen in Mallozzi et al. 
1995). Finally, we use the generated Ep^ak to apply an appropriate bolometric correction 
based on the band function with (a = —1.0 and /9 = —2.0). This bolometric correction ranges 
from a factor of ~3.5 to ~7.1. The resuh is an illustration of the intrinsic distribution of 
the population of bursts on the sky. The point of this Figure is that the distribution covers 
a large area without the Amati relation. 
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Fig. 5. — The Nakar and Piran test for 50 bursts from Amati et al. (2006). This data 
came from a variety of different detectors, including Beppo-SAX, Konus,BATSE, HETE, 
and Swift. All bursts in this sample have known associated spectroscopic redshifts. Of the 
50 bursts, 34% fail the test. This is within the expected failure rate of the Amati relation. 
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Fig. 6. — The Nakar and Piran test for 27 bursts from Schaefer (2007). This data came from 
a variety of different detectors, including Beppo-SAX, Konus,BATSE, HETE, and Swift. All 
bursts in this ample have known associated spectroscopic redshifts. Of the 27 bursts, 41% fail 
the test. This is in within the expected failure rate of the Amati relation and in agreement 
with previous tests on this data (see Schaefer and CoUazzi 2007). 
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Fig. 7.— 1654 BATSE bursts from the future 5B BATSE catalog (Goldstein et al. 2011). 
The Epeak and Sboio data come from the best fit of either a GPL or Band model, whichever 
was significantly better. In addition, we use only the bursts for which the relative error 
on the measurements are 40% or better. This selection fails at a very high rate, with 93% 
violators. The zone covered by the BATSE sample has a very large coverage area, but still 
only a very few bursts are passing. This is particularly condemning, as it hints that any 
future Amati relation will also fail for the BATSE sample. The dotted line represents an 
illustrative line for the trigger threshold and the dot-dashed line represents an illustrative 
model of the Ep^ak detection threshold (see section 4) . 
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Fig. 8. — HETE data from Sakamoto et al. (2005). In total, 44% of the 36 bursts fail below 
the Amati limit line, which is within the expected failure rate. Likewise, 33% of the 12 bursts 
with associated spectroscopic redshift re violators (which again, is within the expected failure 
rate). The bursts with redshift do not appear to be significantly different from those without. 
The HETE data is unique in that it seems to most resemble the original data from Amati, 
even though the scatter around the Amati limit line is very large. The dotted line represents 
an illustrative line for the trigger threshold and the dot-dashed line represents an illustrative 
model of the Ep^ak detection threshold (see section 4). Filled diamonds represent bursts for 
which there is a measured redshift, unfilled circles are bursts for which there is no redshift. 
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Fig. 9. — Swift data from Butler et al. (2007). In total, there are 71 bursts, 82% of which 
violate the Amati limit. This is far beyond the expected value, and thus the Amati relation 
fails for the Swift data. The same conclusion is reached when looking just at the bursts with 
known spectroscopic data, with 76% of those 25 bursts being violators. The bursts with known 
redshift are not different from those without known redshift (see Table 1). The dotted hne 
represents an illustrative line for the trigger threshold and the dot-dashed line represents an 
illustrative model of the Ep^ak detection threshold (see section 4) . Filled diamonds represent 
bursts for which there is a measured redshift, unfilled circles are bursts for which there is no 
redshift. 
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Fig. 10. — Suzaku data. We use only the bursts for which the time-integrated -Bpeafe are 
reported. The fluences are also reported in these notices, typically over the 100 keV to 1 
MeV range. For these fluences, we had to convert them to bolometric fluences using the 
provided spectral parameters. Of the 32 bursts we use, 94% are violators of the Amati limit. 
The 7 bursts that have associated spectroscopic redshift have a 100% violator rate. Again, 
the bursts with known redshift are not different from the overall sample. Filled diamonds 
represent bursts for which there is a measured redshift, unflUed circles are bursts for which 
there is no redshift. 
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Fig. 11. — Combined data from Swift and Suzaku. Krimm et al. (2009) took the raw data 
from both detectors and fit the combined spectra to get a better measurement of Epeak for 
a large sample of bursts. We use their EpeakS as found from the Band function. 38 of the 66 
usable bursts do not have spectroscopic redshifts, of which 86% are violators. These bursts 
have an average log of the energy ratio of 9.63±0.85, whereas the Amati relation requires 
that this average must be less than 9.04. 28 of the bursts have spectroscopic redshifts, 86% 
of those busts arc violators, with an average log of the energy ratio of lO.Olil.Ol. Even 
with the broad spectral range provided by combining the Swift and Suzaku data, the Amati 
relation fails. Filled diamonds represent bursts for which there is a measured redshift, unfilled 
circles are bursts for which there is no redshift. 
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Fig. 12. — Konus bursts. This distribution of bursts has a fairly flat bottom corresponding 
to the trigger threshold for the detector. The very broad energy range of Konus allows for 

Epeak,obs values to be measures from 30 kcV to 2000 kcV. The fraction of bursts violating the 
Amati limit is 73% for the 33 bursts with spectroscopic redshifts and 78% for the the 64 bursts 
without redshifts. With bursts extending down to near the Ghirlanda limit line, the observed 
distribution is clearly not that of the Amati relation plus some ordinary measurement errors. 
In other words, the Amati relation fails for this sample. Filled diamonds represent bursts for 
which there is a measured redshift, unfilled circles are bursts for which there is no redshift. 
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Fig. 13. — Beppo-SAX bursts. These bursts are taken from the Guidorzi (2010) bright 
Beppo-SAX. burst catalog, of which we could use 119 bursts. We find that the fraction of 
bursts violating the Amati relation is 85% for bursts without spectroscopic redshifts, and 
90% for bursts with redshifts. Because these are only the bright bursts, we cannot make 
any commentary as to the distribution of bursts like we do with other detectors. We note 
that despite these being the brightest bursts, this sample has a high violator rate. Since the 
brightest bursts are the most likely to pass the Nakar and Piran test, we can say with some 
confidence that the Amati relation fails for Beppo-SAX bursts. Filled diamonds represent 
bursts for which there is a measured redshift, unfilled circles are bursts for which there is no 
redshift. 
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Fig. 14. — The selection effects. The two selection effects based on the intrinsic distribution 
of the burst population are combined and displayed as contours of burst density. These 
appear as two roughly concave-down parabolas, with each representing a different density 
level. The outside region is shaded darkly so as to indicate that bursts in those regions are 
rare, while the middle region is shaded a light gray to indicate that bursts in those areas of 
the diagram are less common than those in the central area. Of the two detector selection 
effects, the more restrictive is the requirement that the burst be bright enough to measure 
Epeak,obs- We show vcrsious of these detector effects to illustrate how the Amati relation 
bursts could be seen. The lower line illustrates a poor detector threshold (with shading 
below to indicate that no bursts in that area can be measured and placed onto the plot). 
The other line illustrates the result of a detector with both a poor detector threshold and low 
energy range (with shading below it). For a poor detector, the bursts that can be published 
and placed on this diagram are all in the unshaded and unhatched regions. The point of 
this diagram is that the selection effects will force the plotted bursts to roughly he along 
the Amati limit line, and these bursts will then appear to obey the Amati relation. Thus, 
simple selection effects create the Amati relation, at least for some samples of bursts. 



